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Confirmation of strong magnetic field amplification and nuclear 

cosmic ray acceleration in SN 1006 
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Abstract. It is shown that the nonlinear kinetic model of cosmic ray (CR) acceleration in supernova remnants (SNRs) fits 
the shell-type nonthermal X-ray morphology, obtained in Chandra observations, in a satisfactory way. The set of empirical 
parameters is the same which reproduces the dynamical properties of the SNR and the spectral characteristics of the emission 
produced by CRs. The extremely small spatial scales of the observed X-ray distribution are due to the large effective magnetic 
field S d ~ 100 pG in the interior, which is also required to give a good fit for the spatially integrated radio and X-ray synchrotron 
spectra. The only reasonably thinkable condition for the production of such a large effective field strength is an efficiently 
accelerated nuclear CR component. Therefore the Chandra data confirm the inference that SN 1006 indeed accelerates nuclear 
CRs with the high efficiency required for SNRs to be considered as the main Galactic CR sources. 
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A simple-minded calculation (McKenzie & Volk l 19821 of 
the mean square magnetic field (MF) fluctuation (5B) 2 



(6B/B) 2 = M a P c /(p g V s 2 ) 



(1) 



excited by the cosmic ray (CR) streaming instability in a strong 
shock with M a » 1 shows that for efficient acceleration, P c ~ 
p g V 2 , the excitation of resonant Alfven waves is not only an 
important ingredient for the diffusive acceleration mechanism 
as such (Bell [T378l Blandford & Ostriker [T378l The resultant 
(5B/B) 2 » 1 could also imply a considerable obstacle for the 
theoretical description of the process up to the production of an 
enhanced effective MF in the acceleration region (Volk 1984 1. 
In eq^A^a = Vs/c a is the Alfvenic Mach number, where c a 
denotes the Alfven speed, P c is the CR pressure at the shock 
front, and p g V 2 is the ram pressure of the incoming plasma 
flow. The problem of eq^was taken up again recently by Bell 
& Lucek ( 120011 1. referring also to (Lucek & Bell 20(3 and at- 
tempting a nonlinear description of the MF evolution in a sim- 
ple model, with the conclusion that a considerable amplifica- 
tion to what we shall call an effective MF should indeed occur. 
At the same time the selfconsistent treatment of the nonlinear 
time-dependent acceleration equations, both in theoretical and 
computational terms (e.g. Berezhko et al. 1996; Berezhko & 
V61k ll99"71 l. had progressed sufficiently to make quantitative 
modeling of CR acceleration in specific supernova remnants 



SNRs like the historical object SN 1006 possible (Berezhko 
et al. 2002 1, including the so-called injection process (Volk et 
al. 2003 and references therein). This allowed a detailed com- 
parison with radio and X-ray observations which had become 
available in the nineties (Koyama et al. 1995 , Allen et al. 2001 1, 
and allowed in particular a theoretical study of the morphology. 
In fact there are very recent measurements on details of the X- 



ray morphology of SN 1006 with Chandra (Long et al. 2003; 
Bamba et al. 2003 1. We shall show in this Letter that they con- 
firm the amplification of the MF and - for lack of any acceptable 
alternative - its cause, the efficient acceleration of nuclear CRs. 

Nonthermal X-ray observations indicate that at least CR 
electrons are accelerated in SNRs. In SN 1006 there is evi- 
dence that electrons reach energies of about 100 TeV (Koyama 
et al. 119951 Allen et al. 12001b . Also TeV y-ray emission from 
this source has been reported (Tanimori et al. 1998). However, 
depending on the assumed values for poorly known physical 
parameters in the theory (mainly the value of the MF and the 
nucleon injection rate, but also the ambient gas density), the 
observed high-energy y-ray emission of SN 1006 could in prin- 
ciple be predominantly either simply the inverse Compton radi- 
ation of the synchrotron electrons scattering on the microwave 
background radiation as assumed by Tanimori et al. ( 199$}, or 
7T°- decay emission due to hadronic collisions of CRs with 
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gas nuclei (Berezhko et al. 2002). Comparing the synchrotron 
data with the calculated spectrum of the energetic electrons 
Berezhko et al. (2002) have inferred that the existing data re- 
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quire a very efficient acceleration of CR nuclei at the SNR blast 
wave (which converts about 10 % of the initial SNR energy 
content into CR energy) as well as a large interior MF strength 
Bd ~ 120//G. The large value of is expected and possible 
in this case, using the spatially integrated radio and X-ray syn- 
chrotron spectra to empirically determine the MF strength and 
the injection rate in the theory. MF amplification was also one 
of the possibilities discussed later by Long et al. (j2003 1. The 
nonthermal energy is distributed between energetic electrons 
and protons in the proportion 1.5 x 1CT 3 , similar to that of the 
Galactic CRs. From the point of view of acceleration theory 
this is the physically most plausible solution. It can also explain 
the peculiar dipolar structure of the synchrotron emission, with 
the dipole axis parallel to the ambient interstellar MF, as an ef- 
fect of the homogeneity of the large-scale ambient MF (Volk 
et al. 2003i, generally believed to exist in the rather clean inter- 
stellar environment of SN 1006 above the Galactic disk (Long 
et al. 2003). This morphology is fairly common among SNRs 
(Winkler & Long 1997). There is no alternative process with- 
out ad hoc-assumptions in the literature, or a new one which 
we could reasonably imagine, that would amplify the MF in a 
collisionless shock without particle acceleration. Nevertheless, 
the previous observations of the large-scale nonthermal emis- 
sion in the case of SN 1006 did not strongly exclude a scenario 
which we called the "inefficient model", with low interior MF 
Z?d ~ 10 pG in which nuclear CRs would play no important role 
and practically all nonthermal emissions would be of leptonic 
origin (Berezhko et al. 2002 1. As a corollary we shall show here 
that such a low-field scenario is not compatible with the local 
X-ray morphology. 

For the comparison with the new Chandra data we shall 
not present the details of our model. They have been already 
described in the above paper (Berezhko et al. |2"U02). We shall 
rather use some simple analytical approximations, needed to 
qualitatively interpret the spatial distribution of the various CR 
components in SNRs. For the quantitative comparison we shall 
however use the exact numerical results. With some additions, 
the distribution of CRs produced by the spherically expanding 
shock of radius R s and speed V s can be roughly described by 
the steady state 1 -dimensional transport equation for the CR 
distribution function f(r, p, t) 

d 2 f 



df p du df 
dx 2 dx 3 dx dp 



0, 



(2) 



where x = R s - r, and u — V s - w is the speed of the scattering 
medium relative with respect to the shock front, w is the speed 
in the frame of the progenitor star, K(p) is the CR diffusion coef- 
ficient, p is particle momentum, L is a loss term, which we take 
in the simple form L = f/r, where t is the loss time. Within 
this approach we neglect any effects of the shock modification 
due to the CR backreaction. Therefore the velocity profile has 
the form: u(x < 0) = u\ — V s ; u(x > 0) = U2 = U\jcr, where 
<x = 4 is the shock compression ratio. The solution of this trans- 
port equation is /; = /n(p) exp(— \x\/li) (Volket al. 119811 ). where 
/()(/?) = /(* = 0, p) is the CR distribution function at the shock 
front and the scale 

2/C; r » 1-1 



describes the spatial CR distribution in the upstream (i = 1) 
and downstream (i = 2) regions, with simple limits for strong 
(t,- <z Kj/u 2 ) and weak (t,- » ki/u 2 ) losses. The CR distribution 
function at the shock front is determined by the expression 



fa =Ap 9 exp 



•J Pin 



dp<f>(p)/p 



(4) 



l_(_iy^l + 4^. /(M 2 T . ) 



(3) 



where q — 3mi/(mi - M2); <f> — q[K\/(u 2 Ti) + K2/(u\U2T2)]. 

The losses produce two effects. First of all, the universal 
power law spectrum of accelerated CRs / oc p~i has an expo- 
nential cutoff, where the maximum CR momentum p m is deter- 
mined through the condition <p{p m ) - 1 . The losses lead also to 
a reduction of the spatial scales /,■ that is essential in the cutoff 
region p ~ p m , cf. eq@] 

In the collisionless cosmic plasma, losses of nuclear CRs 
due to their interactions with the gas particles or with the am- 
bient fields are negligibly small. The only important effects 
which restrict the proton acceleration are adiabatic cooling in 
the downstream region and a geometrical factor, that is the fi- 
nite shock size. Formally compared with the simplified one- 
dimensional eq|5] the spherically symmetric transport equation 
contains the additional term (2i<i/r)(df/dr) which contains the 
information about the shock size R s . Approximately we can 
write df/dr w -f/h. Therefore this additional term in the 
plane wave approach can be expressed as a loss term L with 
the loss time t\ - R s /(2ui). In the downstream region inside 
the SNR there is in addition also the term (Vup/3)(df/dp) 
which describes CR adiabatic cooling. It can be similarly es- 
timated as a loss term with loss time to ~ R s /V s . As a re- 
sult, the upper proton momentum is determined by the relation 
Kl ( Pm ) = R S VJA, where A ~ 10 (Berezhko [T996l . 

In addition to the above mechanical effects CR elec- 
trons suffer synchrotron losses with a time scale t = 
9m 2 c 2 /(4r^B 2 p), where m e is the electron mass, ro denotes 
the classical electron radius, and B is the MF strength. This 
is the dominant loss effect for our discussion. For sufficiently 
low MF, when the loss time t exceeds the age t of the sys- 
tem, synchrotron losses are not important and the electron and 
proton spectra have exactly the same shape. In this case the 
highest energy electrons with p ~ p m have quite a wide spa- 
tial distribution with scales ~ 0.1R S . For high MF values 
the electron spectrum is restricted to significantly lower mo- 
menta, p e m <K p m , than are protons. In this case accelerated 
electrons occupy only a very thin region around the shock front 
since l\o 0.1/? s for all electron energies, for which t <k t. 
The radio-electrons in the GeV range have t » t and there- 
fore I2 ~ 0.1R S , consistent with observation (Long et al. 2003 1. 
The same is true for 100 TeV electrons which produce X-ray 
emission in the low-field case. 

During the free expansion phase, preceding the early Sedov 
phase in which SN 1006 is at the present time, a Rayleigh- 
Taylor instability of the interface between the shocked inter- 
stellar gas and the ejected gas can occur (Gull 119731 . It can 
make the downstream gas velocity profile w(r) narrower than 
the classical Sedov profile, leading to an increase of the adia- 
batic expansion compared with the above estimate. As a result 
the radio emission profile may become somewhat narrower, 
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Fig. 1. Radial dependence of the X-ray emissivity at different 
X-ray energies. Thick and thin lines correspond to the efficient 
and the so-called inefficient model respectively 



whereas the X-ray emitting electrons are not affected, being 
dominated by the strong synchrotron losses around the shock 
This instability can also locally amplify the MR However, the 
unstable region moves to the deep interior (r > Q).\R S ) as the 
Sedov phase sets in, and dies out. Therefore we do not consider 
these processes any further. In particular they have no influence 
on the CR acceleration which occurs near the shock. 

Our nonlinear model (Berezhko et al. 119961 Berezhko & 
Volk ll997l is based on a fully time-dependent solution of the 
CR transport equation together with the gas dynamic equations 
in spherical symmetry. It yields at any given distance r the mo- 
mentum distribution f(r, p, t) of CRs produced during the SNR 
evolution up to time t. Because of the efficient acceleration of 
the nuclear CR component (mainly protons), consistent with 
the expected injection rate of suprathermal particles into the 
acceleration process (Volk et al. 2003 1, the CRs significantly 
modify the shock structure by their nonlinear backreaction and 
amplify the MF (Berezhko et al. 2002 1. This MF amplification 
on all scales drives the CR scattering mean free path towards 
the particle gyroradius pe, which means isotropic spatial diffu- 
sion. Since the particle distribution is also isotropic this implies 
the Bohm limit in the effective mean MF (Volk 1984): k = 
where v denotes the particle velocity. Fig[2 shows the 
calculated radial dependence of the X-ray synchrotron emissiv- 
ity q v (e v , r) (e.g. Berezinskii et al. 1990 1 for four X-ray energies 
e v = 1, 2, 5 and 10 keV. One can see that the emissivity q v {r) 
peaks at the shock position r = R s because the radiating elec- 
trons have a sharp peak at r — R s due to the synchrotron losses. 
Since within the upstream and downstream emission regions 
the respective MF is approximately uniform, the spatial behav- 
ior of the emissivity q Y (r) is determined by the spatial depen- 
dence of the CR distributions. According to eq|3j the expected 
upstream scale for the energy e y — 5 keV is l\ = 0.01/? s , taking 
into account that the upstream MF is B\ — ct p Bq, the upstream 
medium speed is u\ = V s /cr p , <x p w 2 is the precursor compres- 
sion ratio, V s = 3200 km/s, R s = 7.3 pc, and Bo = 20 fiG is 
the MF in the shock precursor, i.e. the interstellar MF which is 
significantly amplified by the (nonlinear) CR streaming insta- 



bility. This is in good agreement with the numerical results in 
FigE based on a source distance of 1.8 kpc. 

The downstream scale due to the synchrotron losses, ex- 
pected from relations 0, is h ~ s/ K 2 T 2 ~ 0.0037/? s , taking 
into account the values u 2 = V s /cr, B 2 = crZ? , and cr = 6.3 
(Berezhko et al. 2002 1. This agrees with the numerical values 
quite well. Due to the energy-independence of l 2 , the spectral 
form of the downstream q v (e v , r) does not depend on e v (see 
FiglD in agreement with observations (Bamba et al. 2003 1. 

The MF Bi = B(r - R s + 0), just ahead of the gas sub- 
shock at r = R s , and the postshock MF B 2 = B(r = R s - 0) 
are connected by the relation B 2 = cr % B\, where <x s = 3.6 is the 
subshock compression ratio. Since the upstream MF is signifi- 
cantly smaller than the downstream MF, the synchrotron emis- 
sivity undergoes a jump at r — R s , and the emission from the 
downstream region significantly exceeds the upstream emis- 
sion. 

In projection along the line of sight, the radial emissivity 
profile determines the remnant's surface brightness. For the X- 
ray energy interval e\ < e v < e 2 it has the form 

hip) J de v J dxq v (e v , r = ^jp 2 + x 2 ) (5) 

where p is the distance between the center of the remnant and 
the line of sight. It is clear from this expression that, due to 
the shock curvature, the surface brightness profile J v (p) differs 
from the emissivity profile q v (r = p), except for the simple 
case of a plane shock which is parallel to the line of sight. In 
particular, the position p m < R s of the peak value J m = J v (p m ) 
does not coincide with the shock edge p = R s , and the scaling 
values L12 which characterize the spatial brightness behavior 
in the inner (p < p m ) and the outer (p > p m ) regions are not 
simply the downstream and upstream emissivity scales Zi_2. 

The numerically calculated brightness profile for the X-ray 
energy interval between e\ — 2 keV and ei = 10 keV is shown 
in Fig|3] The brightness profile is characterized by the outer 
scale L\ = 0.002/? s = 0.015 pc which comes from the emission 
of the downstream region alone, and from the shock curvature. 
The inner brightness scale is L 2 ~ 1L\ =0.1 pc. 

The sharpest experimental X-ray brightness profiles ob- 
tained by the Chandra observers (Long et al. 2003; Bamba et 
al. 2003 1 are shown in Fig|2 Since the absolute values of the 
measurements are not known, all theoretical and experimental 
profiles are normalized to their peak values. One can see that 
the experimental values agree very well with our calculations. 
All other observed brightness profiles are significantly wider: 
on average L\ = 0.04 pc and L 2 = 0.2 pc (Bamba et al. 2003 1. 
There are several reasons for a broadening of the observed pro- 
file, compared with the ideal case of a spherical shell. First of 
all, it is clear from the observations that the actual shock front 
deviates from a spherical form. The wavy shape of the shock 
front can be produced as the result of a small scale density inho- 
mogeneity of the ambient interstellar medium. Any small scale 
distortion of the spherical emission shell leads to a broadening 
of the observed brightness profile. 

It is important to note that there is direct experimental ev- 
idence that not only the inner part (p < p m ) of the brightness 
distribution J v (p) but also the outer part (p > p m ) is due to 
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Fig. 2. Projected radial dependence of the X-ray brightness in 
the 1.2 to 2 keV (solid) and 2 to 10 keV (dashed) X-ray energy 
interval. Thick and thin lines correspond to the efficient and the 
so-called inefficient model respectively. The Chandra data, cor- 
responding to the sharpest profile, are shown by the histogram 
(Long et al. 2003 1 and the vertical dashes (Bamba et al. l2003> . 



emission from the downstream region r < R s , in contrast to the 
arguments of Bamba et al. ( 2003 1. At the gas subshock the den- 
sity increases by a factor <x s = 3.6 and the gas temperature in- 
creases by a factor of 9.3. Due to this fact the thermal emission 
from the upstream region is small compared to the downstream 
thermal emission. Therefore, if the observed nonthermal X-ray 
emission J v {p) at p > p m is due to the upstream electrons, then 
one should expect that the thermal X-ray radiation has a peak 
value at a smaller distance than p m and drops to almost zero at 
p = p m , if as before p m represents the peak position of the non- 
thermal X-ray emission (in other words the peak of the thermal 
X-rays is expected in this case at smaller distances). However, 
the observed peak positions of the thermal and nonthermal X- 
ray emissions and their shapes in the outer region p > p m are al- 
most identical (Bamba et al. 2003 1. This is a strict confirmation 
that the entire observed X-ray emission comes from the down- 
stream region, exactly as predicted by Berezhko et al. (2002 1. 

In Figs. [0 and |2 we also present calculations for the so- 
called inefficient model. Of course, the number of acceler- 
ated electrons was chosen to be consistent with the observed 
synchrotron emission for the assumed downstream MF B& = 
16 /iG (see Berezhko et al. 2002 for details). Since in this case 
the MF is much lower, so that synchrotron losses play no role, 
the upstream (negligible brightness) and downstream electron 
length scales ~ 0.1/? s are so large (see Fig^ that the bright- 
ness profile is almost two orders of magnitude wider than ob- 
served. The clear conclusion is that inefficient scenarios for CR 
injection/acceleration in SN 1006 should be rejected since they 
strongly contradict the X-ray measurements. 

Note that in order to form such a sharp decrease of their 
distribution towards the center of the remnant one needs a pow- 
erful loss process in the downstream region, independently of 
what specific acceleration process produces the energetic elec- 
trons near the SN shock. Adiabatic cooling, as it was demon- 
strated, is not strong enough. In the diluted ionized plasma the 
only known strong electron loss process is synchrotron losses. 



But this process can influence the 100 TeV electron distribu- 
tion in a significant form only if the downstream MF is as 
high as Z?d ~ 100 jjG, much larger than the typical interstel- 
lar MF. Such a high MF can be created due to CR backre- 
action upstream of the shock. For such a large downstream 
MF B& ~ 100 fiG the number of relativistic electrons con- 
sistent with the observed radio flux from SN 1006 is so low 
(Berezhko et al. 2002 1 that they are not able to create Alfven 
waves with sufficiently high amplitudes SB ~ B. Indeed, tak- 
ing into account that a typical value of the Alfven speed is 
c a = 20 km/s and that for fi<j ~ 100 /iG the pressure of the CR 
electron component f c i « 10~ 3 p g V s 2 (Berezhko et al. 2002 1, we 
have (6B/B) 2 « 0.1 from eq[2 which implies an insignificant 
background MF amplification. The only possibility is efficient 
proton acceleration. In this case their number, consistent with 
all existing data, is so high, P c ~ p g V 2 , that they are able to 
strongly amplify the MF and at the same time to provide effi- 
cient CR scattering approaching the Bohm limit (Bell & Lucek 
2001 1. Such a MF amplification is required by the comparison 
of our selfconsistent model with the synchrotron observations. 
We note that the uniformity of the downstream MF strength 
distribution during the Sedov phase, despite the strong expan- 
sion of the downstream region, is a direct consequence of the 
temporal weakening of the random MF generation at the shock 
Bi °c V s 2 (Bell & Lucek 2001), if we can approximate the MF 
as being isotropic according to the relation B^p^ ^ = const 
(e.g. Chevalier ll974> . 

We conclude that these data confirm the MF amplification 
and the efficient acceleration of nuclear CRs in SN 1006, pre- 
dicted by Berezhko et al. (|2002i. This efficiency is consistent 
with the requirements for the Galactic CR energy budget. 
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